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ABSTRACT 


A  solution  to  Gerard's  differential  equations  for  plastic  buckling  of 
cylindrical  shells  is  found  for  the  case  of  lobar  buckling  under  hydrostatic 
pressure.  An  approximate  formula  based  on  this  solution  is  then  obtained 
for  buckling  in  the  inelastic  region. 

According  to  this  formula,  the  buckling  pressure  is  a  function  of  the 
cylinder  geometry  and  the  secant  and  tangent  moduli  as  determined  from  a 
stress-strain  intensity  diagram  for  the  shell  material.  Agreement  with 
experiments  on  ring-stiffened  cylinders  is  found  to  be  within  4  percent. 

INTRODUCTION 

Experimental  studies  of  the  buckling  of  stiffened  cylinders  under  hydrostatic  pressure 
have  shown  that  collapse  of  the  shell  plating  between  frames  is  frequently  preceded  by  yield¬ 
ing  of  the  shell  material.  This  would  indicate  that  inelastic  shell  buckling  may  be  an  impor¬ 
tant  consideration  in  the  strength  design  of  pressure  vessels,  particularly  when  it  is  realized 
that  residual  welding  and  rolling  stresses  often  induce  inelastic  behavior  at  pressures  well 
below  the  design  strength. 

Inelastic  buckling  of  cylindrical  shells  induced  by  external  hydrostatic  pressure  can 
take  place  in  two  basic  modes:  axisymmetric  buckling,  during  which  circumferential  corruga¬ 
tions  develop  along  the  axis,  and  asymmetric  or  lobar  buckling,  whereby  inward  and  outward 
lobes  appear  alternately  around  the  circumference.  Buckling  of  the  first  type  has  received 
considerable  attention,  but  it  has  usually  been  treated  as  a  failure  due  to  yielding  rather  than 
a  buckling  phenomenon.  Typically,  an  analysis  is  based  on  the  concept  of  an  ideal  material 
which,  at  a  certain  stress  level,  undergoes  an  abrupt  transition  to  the  perfectly  plastic  state. 
The  buckling  pressure  is  determined  not  from  stability  considerations  but  from  the  state-of- 
equilibrium  stresses.  Attempts  to  describe  buckling  of  the  second  type  have  been  rather 
limited  and  have  usually  depended  on  the  intuitive  use  of  a  reduced  modulus  in  place  of 
Young’s  modulus  in  the  elastic  buckling  equation. 

In  recent  years,  however,  advances  in  plasticity  theory  have  made  it  possible  to  ap¬ 
proach  these  problems  more  rigorously.  Investigations  by  Bijlaard,1*2  Ilyushin,9  and 
Stowell,4  among  others,  have  contributed  greatly  to  the  development  of  theory  for  the  inelas¬ 
tic  buckling  of  plates  and  shells.  More  recently,  Gerard5  derived  a  general  set  of  differential 
equations  for  cylinders  from  which  ho  obtained  approximate  solutions  for  torsional  buckling 
and  axisymmetric  buckling  under  axial  compression  for  a  strain-hardening  material. 

Lunchick6  has  since  developed  a  similar,  but  more  exact,  theory  for  axisymmetric  buckling 
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of  ring-stiffened  cylinders  under  hydrostatic  pressure.  From  recent  tests  conducted  at  the 
David  Taylor  Model  Basin,  this  theory  appears  to  be  very  reliable. 

At  the  same  time,  it  appeared  that  the  work  of  Gerard  might  also  provide  a  worthwhile 
approach  to  the  problem  of  asymmetric  or  lobar  buckling  in  ring-stiffened  cylinders,  since  his 
differential  equations  are  sufficiently  general  to  account  for  asymmetric  deformations  under 
hydrostatic  loading.  Work  was  subsequently  initiated  to  obtain  a  solution  to  Gerard's  equa¬ 
tions  for  the  asymmetric  problem.  In  this  report  the  solution  is  derived,  and  an  approximate 
formula  for  inelastic  buckling  is  obtained.  Experimental  data  from  ring-stiffened  cylinders 
are  used  to  evaluate  the  formula. 


ANALYSIS 

PLASTIC  BUCKLING  EQUATIONS 

In  the  Appendix  of  this  report,  the  general  differential  equations  of  Gerard5  for  a  fully 
plastic  cylinder  are  specialized  for  the  case  of  hydrostatic  pressure  loading.  The  buckling 
equations  thereby  obtained  are: 
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where  x  and  a 

u,  v,  and  to 
Es  and  Et 
R 
,h 
v 


are  respectively  the  axial  and  circumferential  coordinates, 
are  the  axial,  tangential,  and  radial  displacements, 
are  the  secant  and’tangent  moduli, 
is  the  radius  to  the  shell  midsurface, 
is  the  shell  thickness, 
is  Poisson’s  ratio, 


D  is  the  bending  rigidity  =  Esh3/ 12(1  -  u2), 

Nx  and  Na  are  forces  per  unit  length  in  the  axial  and  circumferential  directions,  and 
p  is  the  hydrostatic  pressure. 

Since  one  function  recurs  in  all  three  equations,  let 


Et  / 1  du  dv  « 

F  ■  — —  [  —  —  +  —  +  — 

E  s  \  2  dx  da  R 

With  several  differentiations,  u  and  v  can  bo  eliminated  from  Equations  [la]  and  [lb]  so  that 
one  equation  relating  F  and  w  results: 


E,  /  d4F  d4F  d4F  \  d4F  3  Et  d4w 

—  -  +8  -  +  4  -  +  3  -  = - L  - 

Et  '  dx4  dx2da2  ds 4  /  dx4  R  Es  dx4 


With  suitable  differentiations,  Equation  [lc]  can  be  combined  with  Equation  [3]  so  that  a 
single  eight-order  equation  in  tr  is  obtained: 
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where  v4  indicates  the  operator 

A  solution  to  this  equation  can  be  written: 


■  d2  + 

.  dx 2  da2  . 


w  =  A  sin  ka  sin  kx  [5] 

n 

where  k  -  — 

R 


mi t 
L 

L  is  the  U  igth  of  the  shell,  and 
m  and  n  are  integers. 


This  solution  satisfies  the  conditions  of  simple  support  at  the  ends  of  the  cylinder;  i.e.,  that 
d2to 

is?  and - vanish  at  *  ■  0  and  x  -  L.  These  conditions  are  not  unreasonable  for  stiffened 

6x 2 

cylinders  since  it  is  likely  that  the  effective  rotational  restraint  will  be  limited  by  the  forma¬ 
tion  of  plastic  regions  arising  from  high  bending  stresses  near  stiffehers  or  end  supports. 

By  substituting  the  solution  {5]  into  the  differential  Equation  [4],  the  following 
characteristic-value  equation  is  obtained: 


2  C  '  ,R  L 


^  -)-]}- 


[6] 


To  simplify  this  equation  the  following  substitutions  are  made: 

A2  1 


fp 

C 


* 


[7] 


The  equation  is  then  rearranged  so  that  an  expression  for  ppt  the  plastic  buckling  pressure, 
is  obtained: 


2/.0A2  J  1  +  C<f> 

E$  l 

'j  E,hfB4>* 

>+  2 

)  R2  A2 

Stf  [1-0(1-  Z^)] 

[*  * 3  ~T~  ] 

[8] 


•The  subscript  on  /  Indicates  the  atresa  ratio  for  the  plastic  region. 
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ELASTIC  BUCKLING  EQUATIONS 


El  ; 


A  similar  procedure  is  followed  in  obtaining  a  solution  for  the  elastic  case.  When 


is  set  equal  to  unity,  Equation  [4]  reduces  to 


’w  + 


Eh  d4u 
R2  dx 4 


7*  (", 
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dx 2 


+  N. 


6a 2  1 
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which  is  the  Donnell  equation s* 7  for  the  case  of  hydrostatic  pressure  loading.  The  bending 
rigidity  is  now  given  by 


D  = 


Eh3 


12(1  -  v2) 


[10] 


where  uf  is  the  elastic  value  of  Poisson's  ratio.  Substitution  of  solution  [5]  into  Equation 
[9]  yields* 


Pe 


2/e 


r  *2*2  +  _*1" 

Eh  12(1  -  u2)  R2h2 

T+  1-*(1  -fe) 


[11] 


where  fe  is  the  stress  ratio  for  the  elastic  region.  This  equation  could  have  been  obtained 
directly  from  Equation  [8]  by  using  the  elastic  value  for  D. 


•It  is  of  some  Interest  to  compsre  this  result  with  a  similar  equation  obtained  by  Von  Sanden  and  Tolke8  in 
their  comprehenaive  study  of  stability  problems  in  thin  cylindrical  shells.  In  considering  the  elastic  buckling  of 
a  ring-stiffened  cylinder,  they  allow  for  the  variability  of  the  pre-buckling  circumferential  stress  with  the  axial 
coordinate.  Their  buckling  pressure  equation,  in  the  terminology  of  this  report  is: 
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“  midbay  circumferential  stress 


*•  circumferential  stress  at  a  frame 


and  and  Sq  are  determined  from  the  theory  of  Von  Sanden  and  Gunther,^ 


If  <J,  does  not  vary  with  x,  then 


5m  -  50  “  5  “  — 


and  the  buckling  pressure  equstion  reduces  to  Equation  [  1 1], 


1 

2/* 
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MINIMIZATION  OF  EXPRESSIONS  FOR  BUCKLING  PRESSURE 


Equation  [8]  expresses  pp,  the  plastic  buckling  pressure,  in  terms  of  n,  the  number  of 
circumferential  buckling  waves,  and  m,  the  number  of  longitudinal  half  waves.  The  buckling 
pressure  can  now  be  obtained  in  minimized  form.  It  can  be  seen  that  pp  is  effectively 
minimized  with  respect  to  n/m  by  setting 


The  resulting  equation  defines  the  value  of  <f>p  for  which  pp  is  a  minimum: 


[12] 


m4»r4  Et  4 
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~~  18/p  -  3  -  6*p2  (1 


3C^  2 


[13] 


The  corresponding  equation  for  the  minimized  plastic  buckling  pressure  is: 


8>n2ir2£,/p 


c<t>P  t 
—  (3  +  *p  + 


1  + 


3^p2N_, 


l-  3  -  2*p  (1  -  /p)  + 


3^2 


[14] 


which  is  obtained  through  suitable  combination  of  Equations  [8]  and  [13].  Since  pp  in 
Equation  [4]  is  proportional  to  m2,  it  is  clear  that  for  the  minimum  value  of  pp  m  must  be 
equal  to  one  in  all  cases  where  n  is  greater  than  0. 

The  corresponding  minimized  equations  can  be  obtained  for  the  elastic  case.  After 
minimizing  pt  with  respect  to  <f>  in  Equation  [11],  one  obtains  the  equivalents  of  Equations 
[18]  and  [14]  for  the  elastic  case: 


>r4m4  /v/ffA\4  1-2*,  (l-/e) 

12(1- v2)  '  L  /  3-20e(l -/,) 


[15] 
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[16] 


2"2ffl2£/e  /  h\  2/v^A\2r  1 
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These  equations  could  also  have  been  obtained  directly  from  Equations  [13]  and  [14].  Again 
it  is  seen  that  m  must  be  equal  to  1.  For  the  case  where  /  is  equal  to  %  corresponding  to  the 
prebuckling  state  of  stress  in  an  unstiffened  tube,  it  can  be  shown  that  Equations  [15]  and 
[16]  are  exactly  those  given  by  Windenburg  and  Trilling  (Equations  [20]  and  [21]  of  Reference 
10)  in  expressing  the  Von  Mises  buckling  pressure  in  minimized  form. 

Although  Equations  [14]  and  [16]  are  relatively  simple  in  form,  they  contain  the  func¬ 
tion  which  is  not  readily  determinable  from  Equations  [13]  and  [15],  However,  it  is  possi¬ 
ble  to  obtain  an  approximate  expression  for  0  from  a  graphical  representation  of  these 
equations.  Figure  1  shows  plots  of  </>  versus  \J~WK/L  (with  m  equal  to  one)  for  the  following 
cases: 


Elastic: 


Plastic:  — 


1 

2 


0.3 


1 

2 

1 


The  value  of  %  for  Et/Es  was  chosen  as  a  typical  case  for  the  plastic  region.  The  two  stress 
ratios,  and  1,  are  extreme  values  which  should  bound  all  cases  of  practical  interest.  The 
curves  terminate  at  the  line  where  0  is  equal  to  1,  since  this  is  the  case  of  axisymmetric 
(n  =  0)  buckling  for  which  the  minimized  pressure  expressions  no  longer  have  meaning.  These 
curves  suggest  that  a  simple  linear  relationship  between  <£  and  \JlTh/L  might  serve  as  an 
adequate  approximation  for  all  cases.  Use  of  such  an  approximation  implies  that  the  number 
of  circumferential  lobes  is  independent  of  the  material  properties.  After  some  investigation, 
the  equation 


[17] 


represented  by  the  dotted  line  in  Figure  1,  was  chosen  as  a  reasonably  good  approximation.* 
It  is  seen  that  this  line  falls  roughly  midway  between  the  extremes  of  the  curves  presented. 
Use  of  Equation  [17]  in  conjunction  with  Equations  [14]  and  [16]  then  provides  approximate 
expressions  for  pp  and  pf.  It  is  also  helpful  to  make  an  additional  approximation  which 


'Actually,  tha  selection  of  the  factor  1,23  was  somewhat  arbitrary  since  the  buckling  pressure  is  relatively 
Insensitive  to  this  parameter  so  long  as  it  falls  between  1. 0  wnd  1.5.  In  following  an  equivalent  procedure 
for  minimizing  the  Von  Mises  buckling  pressure,  Windenburg  and  Trilling*®  obtained  the  value  1.265. 
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^igure  1  -  Buckling-Mode  Parameter  <f>  as  a  Function  of  \JRh/L 

<f>  has  the  limiting  value  1  corresponding  to  n  =  0  which  defines  axisymmetric  buckling. 


simplifies  Equation  [14]  for  the  plastic  buckling  pressure.  That  equation  can  be  rearranged 
so  that 


3  -  20(1 -/p) 


(1-g) 


[18] 
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where 


C*3(WP)  _ 

2  [3-2^(l -fp)+  -^1] 


[19] 


The  subscript  on  <£p  has  been  dropped,  since  the  single  function  0  is  to  be  used  in  both  the 
elastic  and  plastic  regions.  It  can  be  seen  that  a  will  take  on  its  maximum  value  when  <p  is 
equal  to  1.  For  the  case  E(/Es  =  [a,  a  is  0.091  when  /p  is  y3,  and  0  when  /p  is  1.  Thus  for 
all  cases  where  Et/E s  £  !4>  0  <  a  <  0.091.  Since  Et/Es  will  seldom  be  much  smaller  than  us 
whereas  <£  will  always  be  less  than  1,  the  approximation  that  a  can  be  neglected  will  intro¬ 
duce  only  small  errors.  Equation  [18]  is  thus  reduced  to 


Pp 


UirhY* 

H  \r)  \L  ) 


3  C<t> 


l_8-2*(l -fp)  J 


[20] 


To  examine  the  accuracy  of  the  equations  thus  obtained,  the  results  of  the  approximate  equa- 

P  /  R  \  ^ 

tions  are  compared  with  those  of  the  exact  equations  in  Figure  2,  where  —  [  —  )  is  plotted 

^1'  ^  / 

as  a  function  of  yJun/L ,  p  being  the  theoretical  buckling  pressure.  The  solid  curves  represent 
the  exact  Equations  [14]  and  [16]  for  the  plastic  and  elastic  cases,  respectively,  with  the 
corresponding  values  of  determined  from  the  exact  curves  of  Figure  1.  The  corresponding 
approximate  results,  indicated  in  Figure  2  by  the  dotted  curves,  are  obtained  from  Equations 
[20]  and  [16]  using  the  approximate  expression  [17]  for  <f>.  It  is  seen  that  the  approximate 
method  of  calculation  agrees  quite  closely  with  the  exact  method,  even  though  Figure  1  shows 
wide  divergence  between  the  approximate  and  exact  values  for  0,  particularly  in  the  upper 
range  of  \jlih/L. 


DETERMINATION  OF  SECANT  AND  TANGENT  MODULI 

Before  Equation  [20]  can  be  used,  E s  and  Et  must  be  related  to  the  applied  pressure. 
The  secant  and  tangent  moduli 


[21] 


are  defined  in  the  Appendix  and  are  shown  graphically  in  Figure  3.  For  hydrostatic  pressure 
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Figure  2  -  Comparison  of  Approximate  and  Exact  Methods 
for  Determining  Minimum  Buckling  Pressure 


loading,  and  with  Poisson's  ratio  equal  to  !?,  the  stress  and  strain  intensities  are 


°i  -  (*x2  +  o}  -°xa$)V> 

VT  v  *  •  *  ,} 


[22] 


The  characteristic  stress-strain  curve  of  the  shell  material  is  first  obtained  from  uni¬ 
axial  compression  tests.  In  this  case  ai  and  «,•  are  identical  with  the  axial  stress  and  strain 
(regardless  of  the  value  of  Poisson’s  ratio).  Hence  Es  and  Et  are  readily  determined  from  the 
stress-strain  curve.  In  practice,  it  is  convenient  to  determine  Et  by  drawing  tangents  to  the 
curve. 

Mention  should  be  made  of  one  difficulty  which  may  be  encountered  in  the  interpreta¬ 
tion  of  the  stress-strain  data.  Conventional  strain-measuring  equipment  such  as  an  automatic 
recording  extensometer,  although  adequate  for  measuring  yield  strength,  may  not  be  suffi¬ 
ciently  accurate  for  the  determination  of  Young's  modulus.  Unless  a  high-precision  device 
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Figure  3  —  Typical  Stress-Strain  Diagram 


is  employed,  it  is  best  to  obtain  only  the  relative  values  Ef/E  and  Et/E  from  the  stress-strain 
data  and  assume  a  standard  value  for  E. 

Having  determined  E a/E  and  Et/E  as  functions  of  a,,  one  must  then  apply  them  to  the 
hydrostatically  loaded  cylindrical  shell.  According  to  a  fundamental  hypothesis  of  plasticity 
theory,  the  stress  and  strain  intensities  are  uniquely  defined.4  Thus  by  expressing  hydrostatic 
pressure  in  terms  of  the  stress  intensity,  a  relationship  between  E s,  Et,  and  pressure  will  be 
established.  Since  equilibrium  requires  that  ox  be  equal  to  pR/2h,  only  as  in  Equation  [20] 
remains  to  be  determined.  As  discussed  in  the  Appendix,  os  is  actually  a  continuously  vary¬ 
ing  function  of  x,  whereas  in  this  theory  o s  is  treated  as  a  constant.  Thus  a  single  value  of 
<7S  must  be  chosen,  and  it  is  taken  to  be  the  stress  occurring  midway  between  frames.  Since 
this  is,  generally,  the  maximum  membrane  stress  in  the  shell,  it  might  be  regarded  as  a 
conservative  choice.  However,  it  should  be  noted  that  for  a  material  exhibiting  a  plateau- 
type  stress-strain  curve,  any  other  choice  would  probably  overestimate  the  strength  of  the 
shell. 

In  calculating  o$  it  is  particularly  useful  to  make  the  further  simplifying  assumption 
that  <7S  is  proportional  to  the  applied  pressure  in  both  the  elastic  and  plastic  regions.  This 
assumption  is  reasonable  provided  the  deflections  of  the  shell  remain  small  compared  with 
its  thickness.  Then  o s  can  be  determined  completely  from  the  theory  of  Von  Sanden  and 
Gunther9  with  /  equal  to 


11 


where 


f  sinh  Op  +  sin  6p  \ 
\  cosh  Sp  -  cos  dp  ) 


[23] 


[24] 


Aj  is  the  cross-sectional  area  of  the  frame,  and  b  is  the  faying  width  of  the  frame.  The  sub¬ 
script  p  indicates  that  all  functions  are  given  for  Poisson's  ratio  equal  to  % 

It  will  be  observed  that  the  “beam-column"  effect,  demonstrated  theoretically  by 
Salerno  and  Pulos,*1  is  ignored  in  the  assumption  of  proportional  loading.  This  effect  causes 
a  departure  from  proportional  loading  in  the  elastic  region.  However,  this  departure  is  ordi¬ 
narily  small  and,  in  view  of  the  approximations  already  made,  to  account  for  it  would  be  an 
unnecessary  refinement.  In  those  cases  where  the  effect  is  large  it  can  easily  be  included  in 
the  value  assumed  for  o(>* 

The  stress  ratio  for  Poisson's  ratio  equal  to  %  is  then  given  by 


[25] 


'Additional  daparturas  from  proportional  loading  are  exhibited  by  cylindere  whoae  generator*  are  not  initially 
atralght  Thla  effect  can  be  computed  from  an  analyaia  by  Lunchick  and  Short  ^ 


12 


Solving  Equation  [22]  for  the  applied  pressure,  one  obtains 


*°thfp 


[28] 


A  plot  of  p  versus  o,  from  this  equation  is  a  straight  line  for  the  case  of  proportional  loading 
but  becomes  a  curve  if  the  aforementioned  nonlinear  effects  are  included. 


INELASTIC  BUCKLING 

Although  Equations  [16],  [20],  and  [26]  define  the  buckling  pressure  for  the  elastic 
and  the  fully  plastic  regions,  no  solution  is  given  for  the  inelastic  region  which  lies  between 
these  two  limiting  cases.  However,  by  employing  an  empirical  correction  factor  wherein 
Poisson's  ratio  is  regarded  as  a  variable,  one  can  arrive  at  an  expression  which  reduces  to 
the  proper  limiting  values.  Gerard  and  Wildhorn13  have  found  that  v  can  be  accurately  ex* 
pressed  as  a  function  of  £ a  in  the  inelastic  region  by  the  equation 

-  '  T  -  tG  -  0  1271 


which  reduces  to  %  when  E t/E  is  zero  and  to  ve  when  E  JE  is  one.  Since  Equation  [20]  is 
for  the  fully  plastic  case  where'  v  has  the  value  %,  it  could  be  written 


[28] 


If  %  is  now  replaced  by  v,  a  variable  defined  by  Equation  [27]  and  f  is  replaced  by  /,  a 
function  of  v,  the  equation  for  pc,  the  buckling  pressure  in  the  inelastic  region,  is 


2n2Etf 
30(1  -v2) 


-  W/E, 

1+  ~*\Tt  ~1 


3-20(1-/)  J 


[29] 


With  0  given  by  Equation  [17],  pe  reduces  to  pp  when  £ S/E  is  zero  and  to  pe,  given  by 
Equation  [16],  when  Es/E  is  one.  From  Equation  [29]  it  can  be  seen  that  the  inelastic 
buckling  pressure  depends  on  both  the  tangent  and  secant  moduli. 
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In  determining  the  stress  ratio  /,  at  is  again  taken  to  be  the  stress  midway  between 
frames  as  given  by  the  theory  of  Von  Sanden  and  Gunther,9  but  with  v  a  variable  defined  by 
Equation  [27].  It  is  found,  however,  that  o s  is  practically  insensitive  to  variations  in  v  and 
that  it  is  sufficient  to  treat  /  as  a  constant  which  depends  only  on  tht  Oeometry  of  the  cylinder. 
This  is  the  same  as  assuming  that  of  is  proportional  to  the  applied  pressure.  Since  it  has 
been  found  in  practice  that  variations  in  v  are  small,  ve  can  be  used  for  determining  the  stress 
ratio.  In  this  way  pe  will  still  reduce  to  pe  when  EJE  is  equal  to  one,  and  Equation  [29] 
can  be  written 


where 


and 


2n2Efe  ,k  y  \  L  ) 

w-o'*/ 


».  -  ».  (- 


sinh  0e  +  sin  de 
cosh  df  -  cos  0e 


0t  .  13(1-^],  _ 


[30] 

[16] 

[17] 

[31] 


[32] 


The  subscript  e  designates  functions  based  on  the  elastic  value  of  Poisson’s  ratio. 
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The  inelastic  buckling  pressure  pe  can  now  be  determined  as  follows:  from  the  charac¬ 
teristic  stress-strain  curve  of  the  shell  material,  Ea/E  and  Et/E  are  defined  in  terms  of  oit 
the  stress  intensity.  Hence  pe  can  be  plotted  as  a  function  of  using  Equation  [30]. 
Similarly,  the  applied  pressure  p  can  be  plotted  against  oi  from  Equation  [26].  Since  Equation 
[30]  is  valid  only  when  pe  and  p  are  equal,  the  buckling  pressure  is  obtained  from  the 
intersection  of  the  two  plots. 

Figure  4  illustrates  the  two  general  types  of  material  encountered  in  practice.  The 
first  is  the  strain-hardening  type  (Figure  4a)  which  exhibits  a  continuous  stress-strain  curve. 
Plots  of  Equations  [26]  and  [30]  are  shown  in  Figure  4b,  where  the  buckling  pressure  is 
defined  by  the  intersection  of  the  two  curves.  Figure  4c.  illustrates  the  case  of  an  elastic- 
perfectly-plastic  material.  Since  the  buckling  pressure  abruptly  drops  to  zero  when  the  plastic 
region  is  reached,  collapse  is  simply  defined  at  the  elastic  limit  of  the  material,  as  shown  in 
Figure  4d. 


EVALUATION  OF  THEORY  WITH  EXPERIMENTAL  DATA 
ON  STIFFENED  CYLINDERS 

To  evaluate  the  theory,  results  were  examined  from  previous  tests  of  seven  cylinders 
in  which  asymmetric  (lobar)  shell  failures  were  observed.  The  properties  of  these  cylinders, 
all  of  which  had  external  stiffeners,  are  listed  in  Table  1.  Two  >ere  machined  from  seamless 
steel  tubing  and  had  rectangular  stiffeners.  The  other  five  had  T-stiffeners  and  were  rolled 
and  welded  from  U.S.  Steel  Carilloy  steel  plate.  As  indicated  by  the  symbol  (c)  in  Table  1, 
in  some  cases  specimens  were  taken  from  the  collapsed  cylinder,  whereas  in  others  (unmarked) 
they  were  taken  prior  to  fabrication.  Data  from  uniaxial  compression  tests  of  these  specimens 
were  used  for  the  measurement  of  yield  strength  and  for  the  determination  of  the  secant  and 
tangent  moduli.  Since  all  data  were  obtained  with  an  automatic  recording  extensometer,  which 
was  not  sufficiently  precise  for  the  absolute  determination  of  Young’s  modulus,  only  the  rela¬ 
tive  values  Es/E  and  Et/E  were  determinod  from  these  data.  For  all  cylinders  a  standard 
value  of  30  x  106  psi  was  used  as  Young’s  modulus  E.  Plots  of  oi  versus  «i  for  the  seven 
cylinders  are  shown  in  Figure  5. 

In  Table  2  the  experimental  collapse  pressures  are  compared  with  the  inelastic  buck¬ 
ling  pressures  calculated  from  the  theory  of  this  report  (Equation  [30]).  This  information  is 
also  presented  graphically  in  Figure  6,  where  the  ratio  of  theoretical  pressure  to  the  elastic 
buckling  pressure  pe  (Equation  [16])  is  plotted  against  the  ratio  of  experimental  pressure  to 
Pe- 

Table  2  also  lists  failure  pressures  predicted  by  several  other  criteria.  Elastic  buck¬ 
ling  pressures  were  calculated  from  Equation  [16]  of  this  report,  from  the  theory  o(  Von  Sanden 
and  Tolke,8  and  from  the  theory  of  Von  Mises  using  the  approximate  EM3  formula  [10]. 10 
Shell  failure  pressures  are  given  for  Von  Sanden  and  Gunther  formula  [92a]9  (based  on  simple 
yielding  of  the  exterior  shell  fiber,  at  midbay)  and  for  the  Hencky-Von  Mises  criterion  applied 
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Equation  30 


Elastic  -  Parftctly  Plostic 


Figure  4c  Figure  4d 

Figure  4  -  Graphical  Determination  of  Buckling  Pressure  for  Two  General  Classes  of  Material 
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TABLE  1 


Properties  of  Cylinders 


Shell 

Unsupported 

tm 

Radius 

Length 
of  Shell 

Frame 

Yield 

Cylinder 

Thickness 

■ 

Faying  Width 

Strength 

fl 

h 

L 

Em 

b 

V 

in. 

in. 

in. 

■3-ikhB 

in. 

psi 

T-  2 

38.87 

0.264 

7.24 

1.885 

0.260 

88,000 

(c) 

T-  3 

38.87 

0.260 

8.74 

1.625 

0.260 

108,000 

Welded 

(c) 

with 

T-  6 

26.87 

0.2S6 

7.24 

1.170 

0.260 

115,000 

T-Franes 

T-  2A 

38.87 

0.2S4 

7.24 

0.796 

0.260 

(c) 

103,000 

T-  7A 

.. 

26.87 

0.263 

8.74 

0.683 

0.260 

84,000 

Machined 

0*12 

with 

4.081 

0.0446 

1.06 

0.0248 

0.079 

68,000 

Rectangular 

Frames 

U-22 

4.077 

0.0356 

0.834 

0.0185 

0.068 

70,500 

*A11  yield  strength*  are  defined  et  offset  strain  of  0-002. 


(c)  Specimens  of  shell  material  taken  from  collapsed  cylinder. 

Specimens  for  all  other  cylinders  were  obtained  prior  to  fabrication. 


TABLE  2 

Comparison  of  Theoretical  and  Experimental  Collapse  Pressures 


Cylinder  Number 

Experimental  Collapse 
Pressure 

Inelastic 

Buckling 

Equation  [30]  ( pe ) 

Von  Sanden  and  T6'lke* 

Elastic 

Buckling 

Equation  [16]  (p,) 

Von  Mises10 
(EMB  Formula  [10]) 

Hencky-Von  Mises 
midbay,  midplane 

H 

Von  Sanden  end  Gunther 
Formula  [92A]9 

Welded 

Machined 

T-2 

T-3 

T-6 

T-2A 

T-7A 

U-12 

U-22 

670 

553 

1005 

680 

770 

975 

735 

696 

563 

1016 

705 

748 

938 

734 

930 

631 

1258 

773 

1032 

2014 

1054 

906 

626 

1259 

756 

1010 

1907 

1002 

786 

585 

1180 

705 

995 

1786 

963 

903 

953 

1429 

912 

976 

1081 

777 

695 

742 

1121 

733 

799 

788 

622 
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Compressive  Stress  Intensity,  <r,  in  ksi  Compressive  Stress  Intensity,*,  in  ksi 


0  2000  4000  0  2000  4000  0  2000  4000  6000  8000 

Stroin  Intensity, «.  in  ir/in./in. 

(All  Stroins  Corrected  on  Bosis  of  E*  30  x  I06  psi) 


0  2000  4000  0  2000  4000  0  2000  0  2000  4000 

Stroin  Intensity,  Cj  in  A/in./in. 

(All  Strains  Corrected  on  Bosis  of  E  *  30  *  I06  psi) 


Figure  5  -  Compression  Curves  fot  Cylinder  Materials 

(c)  M (Uriel  taken  from  collapsed  cylinder 

Axle  of  load  cdtreaponda  to  circumferential  direction  in  cylinder 


18 


0.3  04  0.5  0.6  0.7  0.8  0.9  1.0 

Exptrimintol  Pr«*sur«  /pf 

Figure  6  -  Comparison  of  Theoretical  and  Experimental  Collapse  Pressures 


pc  -  Inelastic  Buckling  Pressure  (Equation  [30]) 
pe  =  Elastic  Buckling  Pressure  (Equation  [l6]) 


to  the  membrane  stress  at  midbay,  wherein  failure  is  assumed  to  occur  when  the  Stress  inten¬ 
sity,  oit  computed  from  the  theory  of  Von  Sanden  and  Gunther,  reaches  the  yield  strength  of 
the  shell  material. 

From  Table  2  and  Figure  6,  it  can  be  seen  that  Equation  [30]  is  in  close  agreement 
(within  4  percent)  for  all  cylinders  reported.  In  view  of  the  many  approximations  contained 
in  the  theory,  this  agreement  is  considered  surprisingly  good.  It  is  also  noted  that  for  all 
welded  cylinders  except  T-7A  the  pressures  according  to  Equation  [30]  are  higher  than  those 
observed.  This  is  to  be  expected,  since  such  cylinders  are  weakened  by  residual  stresses 
and  geometrical  imperfections  introduced  during  fabrication,  none  of  which  are  accounted  for 
in  the  theory.  However,  good  correlation  in  these  cases  indicates  that  such  weakening 
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effects  may  not  be  as  severe  as  had  previously  been  suspected.  From  the  uniformity  of  the 
results  there  also  appears  to  be  no  significant  difference  for  the  limited  available  test  results 
between  those  cases  where  test  specimens  were  taken  before  fabrication  and  those  where  they 
were  taken  from  the  collapsed  cylinder. 

CONCLUSIONS 

1.  The  theory  of  this  report,  on  the  basis  of  the  limited  test  data  available,  appears  to 
predict  the  Inelastic  (lobar)  buckling  pressure  with  good  accuracy. 

2.  Final  evaluation  of  the  theory  must  await  additional  experimental  data,  which  should 
include  tests  of  cylinders  with  internal  frames. 
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APPENDIX 


DERIVATION  OF  BUCKLING  EQUATIONS 


In  deriving  general  plasticity  equations  for  cylinders,  Gerard5  defines  stress  and  strain 
intensities  for  Poisson’s  ratio  equal  to  ;£  according  to  the  octahedral  shear  law: 


°i  =  (  °x  + 


2  _ 


o.  + 


3r2)* 


The  secant  and  tangent  moduli  are  then  defined  as 


[33] 


E. 


t . 
i 


dOi 

(It; 


[34] 


With  the  assumption  that  Poisson's  ratio  v  is  equal  to  %  in  both  the  elastic  and  plastic 
region,  the  stress-strain  relations  become 


The  subscript  a  replaces  y  in  Gerard’s  notation  and  refers  to  the  tangential  direction,  as 
shown  in  Figure  7. 

In  treating  the  general  case,  Gerard  considers  a  cylinder  subjected  to  external  loads 
Nx,  Nt,  and  NXI  per  unit  width  and  external  pressure  p.  The  general  differential  equations 
of  equilibrium  are: 


d2u  Al3  d2u  A3  d2u  Al3  d2v  /  Al2  ^3\d2v 
dx2  2  dxda  4  da2  4  dx2  '2  4  'dxda 

A23  d20  dw  A23  dto 

4  da 2  2  R  dx  4  R  da 
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Figure  7  -  Coordinate  System  and 
Cross  Section  of  Stiffened  Shell 
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2  da2  2  dxda  4  dx2  4  dx2  \  2  4  '  dxda 

A j3  d2U  A2  dw  A23  dw 

- +  — - —  ,0 

4  da2  R  da  4  R  dx 


[36] 


d*  to 

D  A. - 

L  dx* 


d*v> 


'13  7T7  +  (^12  +  ^s) 


dx3da 


d*to 

dx2da2 


-  A 


23 


d*w 

dxda3 


+  A. 


lEth  /A  12 

du 

A  33  du 

^23 

dv 

ZR  '  2 

dx 

4  da 

4 

dx 

d2v> 

d3v 

d2V> 

+  Nt -  +  2  Nxt 

dx2 

- +  N, 

dxda 

da2 

+  p 

da  R 


d4«-| 

~dS  J 
") 


In  this  analysis,  long  cylinders  which  buckle  in  the  oval  («  -  2)  mode  will  be  excluded  from 
consideration,  since  the  differential  equations  are  derived  with  certain  approximations  which 
are  valid  only  when  n1»l. 


The  bending  rigidity  for  v  equal  to  %  is 


D  - 


Eth* 


Qerard  defines  the  plasticity  coefficients  as  follows: 


Ai  m  1  -  -4- 


(  a°i 
*2  -1--T 


i43  -  1  -  ar  2 


A2l  m  A12  m  1  2" 


^31  “  ^13  “a<7xT 

A32  “  ^23  “aff«T 


[87] 


[38] 


where 


In  the  case  of  uniform  hydrostatic  pressure,  it  is  readily  seen  that. 


[39] 


#*• 
h  “° 

[40] 

Nx  pR 
*  h  "  2A 


On  the  other  hand,  ot  is  not  so  easily  disposed  of.  In  deriving  his  equations,  Gerard  con* 
aiders  N g  as  independent  of  z,  an  assumption  which  is  correct  in  the  two  cases  (axial 
compression  and  torsion)  for  which  he  obtained  solutions,  but  one  which  is  clearly  not 
correct  in  the  case  of  hydrostatic  pressure  if  the  ends  of  the  cylinder  are  restrained. 
However,  without  such  an  assumption  it  would  be  difficult,  if  not  impossible,  to  solve  the 
problem  by  means  of  differential  equations.  Furthermore,  in  order  to  simplify  the  plasticity 
coefficients,  the  additional  assumption  is  made  that 

N,  pR 


whereby  the  stress  intensity  becomes 


°i  m  [42] 

and  the  resulting  plasticity  coefficients  are 


Although  it  would  logically  follow  that  .Vf  in  Equation  [36]  should  be  replaced  by  pR,  this 
substitution  is  not  necessary  for  the  solution  of  the  differential  equations.  Instead,  by 
retaining  Nt  as  an  arbitrary  factor,  a  more  general  solution  can  be  obtained.  The  load  Nt 
can  later  be  determined  from  appropriate  theory,  as  described  in  the  body  of  this  report.  The 
resulting  differential  equations  are: 

Et  /  i  d2u  d2v  i  dio  \  3  d2u  1  d2u  ^  d2v 

- 1  -  —  +  -  +  —  —  )+ - + - +  -  -  «  0  [la] 

2 Et  \  2  dx2  dxda  R  dx  /  4  dx2  4  da2  4  dxda 


Et  (  j  d2u  d2V  i  dto  \  j  d2v  j  dau 

- (  «  -  +  -  +  -  —  1+7  -  +  -  -  -  0  [lb] 

Et  \  2  dxda  da2  R  da  /  4  dx2  4  dxda 

4Eth  /£,  \  /  j  du  dv  u»\  T  E,  /j  d4 w  d4w 

SR  'E/'2d*  da  r)  Le.U^  dx2ds2 

d4u>\  o  <J4w  d4u>  "I  d2w 

+ - )  +  -  -  + -  +  AL  - 

da4  I  4  dx4  dx2da2  J  dx2 
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